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Abstract: This paper proposes a mathematical justification of the phenomenon 
of extreme congestion at a very limited number of nodes in very large networks. 
It is argued that this phenomenon occurs as a combination of the negative 
curvature property of the network together with minimum length routing. More 
specifically, it is shown that, in a large n-dimensional hyperbolic ball B of radius 
R viewed as a roughly similar model of a Gromov hyperbolic network, the 
proportion of traffic paths transiting through a small ball near the center is 
9(1), whereas, in a Euclidean ball, the same proportion scales as 9 ( J?n 1 _ 1 ) . This 
discrepancy persists for the traffic load, which at the center of the hyperbolic 
ball scales as vol 2 (B), whereas the same traffic load scales as vol 1+ ™ (B) in the 
Euclidean ball. This provides a theoretical justification of the experimental 
exponent discrepancy observed by Narayan and Saniee between traffic loads 
in Gromov-hyperbolic networks from the Rocketfuel data base and synthetic 
Euclidean lattice networks. It is further conjectured that for networks that do 
not enjoy the obvious symmetry of hyperbolic and Euclidean balls, the point of 
maximum traffic is near the center of mass of the network. 



1 



1 Introduction 



One of the most important challenges in large and wide area networks is to 
overcome the traffic congestion problem. The queuing feature at routers can 
create a logic bottleneck between two users. Correspondingly, insufficient band- 
width on the physical links between routers is a contributor to congestion. The 
current congestion control technologies in communication networks are based 
on the feedback from the congested node to the source to slow down the packet 
flow rate, such as bidirectional congestion control and Random Early Detection 
(RED). However, these technologies can only be applied once the congestion 
has happened to some degree, and it is only based on the local point of view of 
some queue overflow along the source-to-target path. 

From the large scale point of view, it has been experimentally observed that, 
on the Internet and other networks, traffic seems to concentrate quite heavily 
on some very small subsets. The major result of this paper is to show that the 
deeper reason behind this congestion in the large scale is the combination of the 
least cost routing together with the negative curvature of the network. 

Roughly speaking, a network is negatively curved if its graph can be approx- 
imated by a negatively curved Riemannian manifold. In practice, however, it 
has become customary to check the Gromov Thin Triangle Condition (TTC), 
meaning that the least cost paths between three vertices "arch inside" the trian- 
gle, giving it a "thin" external appearance. The connection between negatively 
curved surfaces and the TTC can be understood on the basis that a triangle 
drawn on a negatively curved surface has the sum of its internal angles < n, 
giving it a "thin" appearance. The formalization of this equivalence is, however, 
must harder and is known under the Bonk-Schramm theorem [3]. 

Over the past few years, there has been mounting evidence that many (wired 
and wireless) communication networks are negatively curved [18, 11, 13, 12]. Hi- 
erarchical networks have been found to have "hidden hyperbolic structure" [17]. 
In view of the above graph-manifold identification, we can certainly argue on a 
negatively curved n-dimensional Riemannian manifold. To simplify the exposi- 
tion in this Introduction, we invoke the nontrivial fact, proved in the Appendix, 
that congestion in very large hyperbolic ball models does not depend on the 
dimension. Therefore, it is legitimate to argue on a negatively curved surface, 
as we will do in the remainder of this Introduction. 

A well-known feature of least length paths on a negatively curved surface 
is the fact that two geodesies starting from points within an arbitrarily small 
disk eventually diverge exponentially. In this paper, we somehow reverse that 
argument and show that least length paths departing at remote points of a 
convex subset X of the surface "converge" to a single point where the "density 
of geodesies" is the highest. The latter is the point of maximum congestion. We 
conjecture that this congestion point is near the center of mass of the convex 
subset, and prove that the maximum congestion point and the center of mass 
coincide for very large disks embedded in the surface. 

Naturally, there might be some questions as to the relevance of this idealized 
analysis to real networks. In fact, some recent experimental results [20] have 
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confirmed the validity of our analysis. The network congestion metric adopted 
in [20] is the betweennness centrality [24]. The continuous geometry traffic 
congestion metric adopted here is the average length of the geodesies in a small 
convex subset X. We call the later load measure, A t (X), because, if X is 
thought as a sub- network, A t (X) is meant to be the number of packets in that 
sub-network. Consistently with the major theme of this paper, we prove that, 
if X is a disk of fixed radius, A t (A) is maximum when the disk has its center at 
the center of mass. But more importantly, we show that, in a large disk Br(0) 
of hyperbolic surface, maxxc_B H (o),aroa(X)=i ^-t(X) scales as (area(i?^(0))) 2 . On 
the other hand, it was shown in [20] that, for networks for the Rocketfuel data 
base [22], the maximum betweennness centrality scales as TV 2 , where N is the 
number of vertices. Since N can be interpreted as the area, this confirms the 
validity of our theoretical model. 

An outline of the paper follows: In Section 2, we set the basic traffic met- 
rics in graph models. In Section 3, a simple concept of negatively curved planar 
graphs is introduced and the general facts about congestion in negatively curved 
graphs are proposed as conjectures. Immediately thereafter, in Section 4, we 
show that our conjectures hold in selected planar graph examples. In Section 5, 
we propose the more general concept of Gromov-hyperbolic graphs and we in- 
voke the Bonk-Schramm theorem to argue that traffic in Gromov-hyperbolic 
graphs can be analyzed on negatively curved Riemannian manifold models. In 
Section 6, we formulate our major results dealing with traffic in subsets of HP. 
Probably the most significant results are Theorems 3 and 2, proving that traf- 
fic in a small subset X of large hyperbolic and Euclidean balls -Br(O) scale 
as vo1(.Br(0)) 2 and vol(.Bft(0)) 15 , respectively. (The proofs are relegated to 
the Appendix, though.) Theorem 4 and 5 formulate the basic minimum iner- 
tia versus maximum traffic issues, proving the conjectures for hyperbolic balls 
X. Finally, in Section 7, we show that no matter how theoretical our models 
are, they turn out to be surprisingly accurate at predicting asymptotic traffic 
distribution in realistic networks. 

2 Traffic metrics 

Let G = (V, E) be a (connected) graph specified by its vertex set V and its edge 
set E and endowed with a (symmetric) distance function d : G x G — > M. + . A 
path p(s,t) from s to t is a continuous map [0,/] — ► G such that p(s,t)(0) = s 
and p(s,t)(l) = t. The weight of an edge e = xy is defined as w(e) = d(x,y). 
The length of the path is defined as £(p(s, t)) = X) e c P ( s t) w ( e )- A geodesic [s, t] 
is a path such that i ([s, t]) < t (p (s, t)), for all paths p(s, t) joining s to t. 

The traffic on the graph is driven by a demand measure Ad : V x V — > M + , 
where the demand A<j(s, t) is the traffic rate (e.g., number of packets per second) 
to be transmitted from the source s to the destination target t. Assume that 
the routing protocol sends the packets from the source s to the target t along 
the geodesic [s,t] with probability 7r([s,t]). It is indeed customary as a load 
balancing strategy to randomize the Dijkstra algorithm so as to distribute the 
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traffic more evenly [19, 5]. Under this scheme, the geodesic [s,t] inherits a 
traffic rate measure r([s,i]) = Ad(s, t)ir ([s, t\). An edge e laying on the path 
[s,t] inherits from that path a traffic r([s,t]). Aggregating this traffic over all 
source-target pairs and all geodesies traversing the edge e yields the traffic rate 
sustained by the edge e, 

r ( e )= E E r (M)= E E MMMM)- 

(s,t)eVxV [s,t]De {s,t)£VxV [s,t]De 

Observe that r(e)/(e) can be interpreted as the traffic load, that is the number 
of packets in the edge e. 

This paper is essentially concerned with existence of a sub-network that has 
extremely high traffic load. Formally, given a connected subgraph X C G, 
containing some edges, we define its traffic load to be representative of the 
number of packets in it: 

A*(A) = £ e€Jf *(e)r(e) 

= Es.tev (Eee[s,t]nx w ( e )) MMMM) (1) 
= E s ,tevU[s,t]nX)A d (s,t)n([s,t}) 

The above definition does not allow X to be reduced to a vertex, but we 
can identify the smallest A t -mcasurable neighbor of a vertex. Define the star of 
a vertex, star(x), to be the smallest subgraph of G containing x; let /(star(x)) 
be the sum of the lengths of set of edges abutting to x; the latter vertex set is 
denoted as E x — {xy : y G V, xy 6 E}. Then consider 

A t (star(a;) ) E e gi; x A t (e) ^ A t (e) \ ^ 
/(starC*)) E e , E J(e) " /(e) ' T ^ 

The inequality is a well known fact, and equality is achieved for the number of 
hops metric, that is, w(e) = 1, Ve G i?. Since the traffic in any edge connected 
to x must be "serviced" by the "hub" x, the interpretation of t(x) is the traffic 
rate sustained by x. Therefore, the fundamental question addressed by this 
paper can be reformulated, from a very local point of view, as that of existence 
of vertices with very high traffic rate. The above string indicates that such 
vertices with very high traffic rates can be sought via the load measure A t . The 
latter aspect will play a crucial role in Section 6. 

The traffic t(x) can be computed by counting the number of geodesies [s,t] 
having a: as a vertex. One should consider two classes of such geodesies, though: 
those that traverse x, that is, x ^ s, x ^ t, and those that either start at x, 
[x,t], or terminate at x, [s,x]. The first ones have to be serviced twice, once 
in the input queue, once in the output queue; the others have to be serviced 
just once. There are 2(7V — 1) geodesies starting or terminating at x. The 
betweenness centrality (3 C of a vertex v is defined, for a uniquely geodesic 
graph, as the number of geodesies that have v as a vertex (this includes those 
geodesies starting at v or terminating at v). For uniformly defined demand, 
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that is, Ad(x,y) = l,Vx ^ y, the connection between the traffic rate and the 
betweenness is easily seen to be 

t(v) = 2{/3 c (v) - 2{N - 1 )) + 2{N - 1) = 2/3 c («) - 2{N - 1) 



3 Basic conjectures 

For the sake of simplicity, we introduce a network curvature concept restricted 
to planar communication graphs and based on Alcxandrov angles [4]. 

Let (abi — a&deg(a)+i> a ^2, ■ a^dcg(a)) be a cyclic ordering of the set of edges 
attached to the vertex a. A geodesic triangle is defined as Aabc — [a, b] U [6, c] U 
[c, a]. The Alexandrov angle a k at the vertex a of the geodesic triangle Aab k b k+ i 
is 

_ cog _i d(a, b k f + d(a, b k+1 ) 2 - d(b k , b k+1 ) 2 
2d(a,bk)d(a,bk+i) 
Then the (Gauss) curvature at the vertex a is defined as 

n(a) = dc ^ ^ Qfc 

Z)fe=i area(Aa6 fc 6 fc+ i) 

where area(Aafefc6fe + i) denotes the area of the geodesic triangle Aab k bk+i, easily 
computable from the distances via Heron's formula, ft is easily seen that, for 
the number of hops metric (w(e) = f), we have a k = 7r/6; therefore, k(o) < 0, 
k(o) = 0, or k(cl) > depending on whether deg(a) > 6, deg(a) = 6, or 
dcg(a) < 6, respectively 

This simple definition is introduced to construct some easily understood 
illustrative examples. In the main body of the paper, though, we will use the 
Gromov definition of negative curvature for graphs (see Sec. 5.) 

Definition 1 ([14]) The moment of inertia of a connected weighted graph 
G with respect to a vertex v is defined as (v) = J^v^V d 2 {v, v l ). 

Since the edges are treated as massless, this concept refers more to the underly- 
ing metric space structure (V, d) than to the weighted graph structure (V, E, d). 
Observe that this inertia may be infinite. 

Definition 2 ([14]) A center of mass or centroid of the weighted graph G 
is defined as a vertex that achieves the infimum of the inertia: 0^(cm(G)) = 
inf„ e y 4> G {v) < 00. 

If the graph is infinite (|V| = 00) , this definition requires existence of a vertex v 
such that <f>Q (v) < 00. If we relax the minimum to be anywhere on the graph, 
it will in general be achieved on an edge, but restricting it to vertices makes it 
easier to relate it to the traffic. 

We can now pose our conjectures. 



5 



Conjecture 1 (Negative curvature) Consider a large but finite (\V\ < oo) 
negatively curved graph G, subject to uniformly distributed demand. Then 

1. There are a very few nodes v that have very high traffic rate t(v) as mea- 
sured by p c (v); furthermore, the vertices of highest traffic rate are in a 
small neighborhood of the vertices of minimum inertia. 

2. If the graph has a symmetry group that fixes some point 0, this point 
achieves the unique minimum of the inertia and the maximum of the traffic 
rate. 

The first part of this conjecture is illustrated on planar graphs in Section 4. As 
a model of a large negatively curved graph, it is tempting to take the full hy- 
perbolic space HP. The latter is (globally) afhne symmetric [8, Chap. IV], [15, 
Sec. 6.2], [16, Chap. XI, Example 10.2], and the transitivity of the symmetry 
group [16, XI, Th. 1.4] would make the traffic uniformly distributed, if it weren't 
for the lack of convergence of the traffic function for uniformly distributed de- 
mand on infinite space. Restricting the network model to finite subsets of H™ 
breaks enough of the symmetry to create traffic spikes, even asymptotically as 
the size increases to infinity. An estimate of how close the center of mass and 
the point of maximum traffic is seems to be beyond our reach for general net- 
work models, at least at this stage. However, for hyperbolic balls, we prove that 
the two coincide and, as major result, that the sharpest traffic spike scales as 
vol(G) 2 . 

Conjecture 2 (Nonnegative curvature) Consider a large but finite (\V\ < 
oo ) nonnegatively curved graph G subject to uniformly distributed demand. Then 

1. Both the traffic and inertia functions r and f3 c are more evenly distributed 
than in the case of a negatively curved graph. 

2. If the graph has a vertex transitive symmetry group, then both the traffic 
and the inertia are uniform. 

Even though this conjecture asserts that for both zero curvature graphs (e.g., 
Euclidean lattices) and positively curved graphs, the traffic is more smoothly 
distributed than in the case of negatively curved graphs, there is a significant 
difference between the two cases. A Euclidean lattice graph could be infinite, 
while positively curved (cubic) graphs are finite by Higuchi's theorem [9, 23, 
21, 6]. Positively curved graphs need not be truncated and hence enjoy more 
symmetry than truncated Euclidean or hyperbolic graphs. This point is easy to 
illustrate on the 1-skeleta of the boundaries of the Platonic solids, all of which 
are positively curved. (For those Platonic solids that have triangular faces, the 
above curvature formula applies; for the other Platonic solids, the more general 
Higuchi-Mohar-DeVos formula [6] should be applied.) The transitivity of the 
symmetry groups of the Platonic solids, together with the uniform distribution 
of the demand, implies the uniform distribution of the traffic on the boundary 
graphs. For more general positively curved graphs, this conjecture is proved 
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using a Ricmannian manifold model with its curvature bounded as < k\ < 
k(x) < fc|. Here the symmetry is broken by the nonisotropic curvature 

Truncating E™ to secure convergence results in symmetry breaking, creating 
traffic "bumps," not as sharp as those of the negatively curved case, but sharper 
than those in the positively curved spaces. This statement is quantified by 
showing that the maximum traffic scales as vol(G) 1+ ~. 

4 Some examples 

We consider two simple examples: The first one (symmetric graph) highlights 
the worsening of the congestion as the curvature decreases. The second one 
(asymmetric graph) highlights the relationship between the curvature and the 
inertia. 

4.1 Almost symmetric graph: congestion versus curvature 

We construct a graph from a single vertex (#1) followed by the addition of seven 
neighbors (#2-#8) in a counterclockwise sense, as shown in Figure 1, top. We 
then proceed from vertex #2 and add vertices #9, #10, #11, #12, so that 
vertex #2 has valence 7. We then add neighbors to vertex #3, and proceed 
recursively, until we obtain the graph shown in Figure 1, top. With 7 neighbors 
for each vertex, except the boundary ones, the graph is negatively curved, be- 
cause ^2i—i on = 7 x -| > 2ir. The same construction can be done to generate 
graphs of valence 6 (vanishing curvature) and 8 (negative curvature). The plots 
of Figure 1, bottom, clearly demonstrate that the maximum congestion worsens 
as the curvature becomes more negative. 

4.2 Asymmetric graph: congestion versus inertia 

Here we consider a graph of valence 7, constructed the same way as in the 
preceding case, except that some "appendices" have been deliberately added 
to make the graph asymmetric. The results are shown in Figure 2. Because 
the graph is negatively curved, and as conjectured, the traffic has its maximum 
precisely at the point where the inertia is minimum. Observe that the traffic 
and the inertia are "in opposite phase." 

5 From graphs to Riemannian manifolds 

As said in the Introduction, it has become standard routine in complex networks 
to check their Gromov negatively curved property as a way to validate argu- 
ments based on such models as the Poincare disk. For the sake of mathematical 
accuracy, we formulate the basic definition and result. 
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Figure 1: Top: Almost symmetric negatively curved network of valence 7. Bot- 
tom: Traffic function r versus vertex number; clearly, congestion at centroid 
increases as curvature becomes more negative. Also shown is the traffic func- 
tion after curvature-based load balancing (see [f9] for details). 

Definition 3 A geodesic metric space (G,dc), for example a graph, is said to 
be Gromov hyperbolic if there exists a 5 < oo such that every geodesic trian- 
gles Aabc has an inscribed triangle Axyz, x £ [b, c],y 6 [a, c],z £ [a, b], of a 
perimeter not exceeding S, that is, d(x, y) + d(y, z) + d(z, x) < S. 

This definition makes sense only for infinite graphs. Even though there exists a 
Gromov concept for finite graphs [11], we really do not need the latter here, as 
our congestion analysis is asymptotic for very large graphs (or manifolds). 

Theorem 1 (Bonk-Schramm [3]) Let (G,da) be a Gromov hyperbolic met- 
ric geodesic space with bounded growth at some scale. There exist an integer n, 
a convex subset D C H n , constants A, k, and a map f : G — > D such that 

\Xd G (u,v)-d D (f(u)J(v))\<k, Vu,veG 

and sup^g^, do(x, f(G)) < k. 
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Figure 2: Traffic and inertia versus vertex number on a planar graph of valence 
7; the graph has been made asymmetric so that the conjecture could be verified 
on cases where the centroid does not coincide with the obvious symmetry center. 
The centroid occurs at vertex #14. Also observe the heavy traffic at vertex #4, 
because it is the gateway between the top appendix and the right side of the 
network. Observe the "opposite phase" phenomenon between the inertia and 
the traffic, clearly seen for vertex #86. 

This theorem makes precise the somewhat loose statement made in the intro- 
duction regarding graph-manifold identification. Clearly, G, after a scaling A, 
can be identified with D via / and subject to a bounded error k. The condition 
of bounded growth at some scale is satisfied, for example, when G is a finite 
valence graph. 

Since this graph-manifold identification entails a bounded error, large scale 
problems on graphs can be mapped to more manageable continuous geometry 
problems on manifolds. 
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6 Differential geometry proof of conjectures 



To justify our numerical results related to congestion on planar graphs of uni- 
form valence of 7, 8, 9, we could develop a Poincare disk model (see Fig. 3, 
right panel). The latter is a faithful model in the sense that the graphs of va- 
lence 7, 8, 9 are quasi-isometric to the Poincare disk. Recall that the Poincare 
disk D = {z = x + jy € C : \z\ < 1} inherits its hyperbolic structure through the 
metric ds 2 = ^ d * d ^yi ■ The latter leads to the area element darea = 

In the spirit of generalization, we develop models in the hyperbolic space HP, 
even in a Riemannian manifold M with arbitrary Riemannian metric ds 2 and 
volume form dvol, provided its sectional curvature bounded as — k 2 < n(x) < 
— fc 2 < 0. To emphasize the role of the curvature, we also look at the Euclidean 
space E 2 as a model of graphs of valence 6 and further generalize the results to 
E™. 



6.1 Hyperbolic versus Euclidean traffic load 

Consider, in a Riemannian manifold M, a large ball Br(0) of radius R with 
its center at the origin together with a convex subset X C Br(0). The traffic 
load in X is defined as (compare with (1)) 

A*PO= ff l(Xn[s,t])dA d (s,t) 

J J(s,t)eB R (0)xB R (0) 

As for the graph model, in this Riemannian context, the demand is uniformly 
distributed in the sense that dAd(s, t) is the product volume dvol(s)dvol(t) . The 
normalized traffic load in X is defined as follows: 



X t (X) := 

= vol(S„(0)) 2 JJ(s,t)eB R (0)xB R (0) 



vol(B„(0)) 2 



IL,t)eB R (0UB n (0) 1 ( X n t s > *]) dA ^ *) 



We take the hyperbolic plane H 2 as a roughly similar model of an infinite 
negatively curved planar graph [9] and carry over the analysis to n dimensions. 
In order to secure convergence of the traffic load, we restrict ourselves to the 
finite domain Br(0) C HP, representative of a large but finite negatively curved 
graph G, and then we do the asymptotic analysis as R — > oo. 

We now state our two major results. The proofs are in the Appendix: 

Theorem 2 If -Br(O) and X = B r (0), r« JJ, are concentric balls in HP and 
if dkd(s,t) is the product volume measure, we have 

A t (B r (0))xci(n)r n 

where Cfe(n) > 0. Furthermore, in a Riemannian manifold of curvature bounded 
as —k 2 < k(x) < —k\ < 0, \fx G M, the asymptotic traffic load is bounded as 
c kl {n)r n <\t(X)<c k2 (n)r n . 
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Figure 3: Traffic load in Euclidean space and Poincare disk, the latter taken as 
a model of H 2 . The radii of the balls in the complex plane model are written 
f,R and the corresponding hyperbolic measurements are r = tanh _1 (f), R = 
tanh -1 ^). 

Theorem 3 If -Br(O) and X = B r (0), r -C R, are n- dimensional Euclidean 
balls (vanishing curvature) and if dAd(s,t) is the product volume measure, we 
have 

\ t (B r (0)) x co(n)-^ 

for some constant Co(n) > 0, independent ofr,R. 

Probably the most important conclusion to be drawn from the preceding 
two theorems is that the normalized traffic load in the small hyperbolic ball 
B r (0) remains bounded from below as R — > oo, whereas the same normalized 
traffic but in the small Euclidean ball goes to zero as R — > oo. Figure 3 pro- 
vides an intuitive explanation as to why this discrepancy happens. Because the 
hyperbolic geodesies are "arched" towards the center where the small ball lies, 
their average length in the small hyperbolic ball is much larger than in the small 
Euclidean ball. 

Figure 3 gives a clue about elementary proofs of Theorems 3-2 in 2 dimen- 
sions. Clearly, the natural parameterization of the s(x,y), t(x',y') points in 
the traffic load integral is via polar coordinates; the difficulty is to compute the 
Jacobian from dxdydx' dy' , ( 1 _( x 'i + y%fy d ^^ x ,'i +y riys ) -i , resp., to the area squared 
element in polar coordinates of Euclidean, hyperbolic, resp., spaces. The de- 
tails are available in [19]; this elementary but more explicit proof yields a specific 
value for co(2) = 1/tt. The proof in the appendix is much more conceptual. 

Observe that, besides its definition as the traffic load, A t could be interpreted 
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in another way. Instead of (2), consider the following: 



1 



// 

J Je 



Ixn[ s ,t](s,t)dA d (s,t) 



vo\(B R (0)y 



B R (0)xB R (0) 



where /x n | S]f ](s,i) = 1 if X n [s,t] ^ and otherwise. The above is clearly 
the proportion of communication paths transiting through X . If X is between 
two balls, that is, B ri (0) C X C B r2 (0), the above is of the order of X t (X)/ri. 
Thus Xt(X) can be interpreted as the percentage of traffic passing through X. 

6.2 Minimum inertia 

Let D C M be a convex domain of finite volume in a Ricmannian manifold with 
distance d(-, •). For x £ D, define (compare with Def. 1) 



For p = 2, the above is called moment of inertia of D relative to x. Next, 
assume there exists a point cm(D) £ D such that (compare with Def. 2) 



For p = 2, cm(D) is called a center of mass or centroid of D (sec [14, Def. 
3.2.1]). This concept was apparently introduced by Elie Cartan [2, p. 47]. 

Arguments related to existence and uniqueness of the center of mass rely on 
strict convexity of x ^ d(x,y) p and x i— ► 4>^{x). In negatively curved spaces, 
p > 1 suffices, whereas in nonpositively curved spaces (e.g., Euclidean spaces) 
the stricter condition p > 1 is required. If <fi^ is strictly convex, uniqueness is 
guaranteed [14, Lemma 3.1.1]. 

Theorem 4 Let M be a complete Riemannian manifold with its curvature 
bounded as —k% < k(x) < —k\ < 0. Then the inertia of -Br(O) C M rela- 
tive to the point x, 



has a unique minimum; furthermore, for k\ — k 2 , this minimum is x = 0. 

Proof. Obviously, M is a complete Riemannian manifold of nonpositive cur- 
vature, and hence it is a Busemann Non Positively Curved (NPC) space [14, 
page 45]. Define the measure = vo\(-)I Br ^(-), where Ib r (o) denotes the 
indicator of -Br(O) C M. Obviously, the measure /x has finite support and 
a(B R (0)) < oo. Therefore, by [14, Th. 3.2.1], 




4> ip) (cm(£>)) = inf (p) ( x) < oo 
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is finite and has a unique infimum. Furthermore, if k\ = ki = 1, M has 
orientation preserving isometry group SO + (n, 1). Under the subgroup that fixes 
0, .Br(O) is invariant, and therefore (j)^ has the same symmetries. Therefore, 
the only way to secure uniqueness of the infimum is x* = 0. 



6.3 Maximum load 

Theorem 5 Restricted to convex subsets X of -Br(O) c H™ with the same 
hyperbolic area, X t (X) reaches its maximum for a ball centered at the origin of 
B R (0). 

Proof. The proof is a corollary of the proof of the Appendix. First, observe 
that /^^(x) is spherically invariant, that is, /Ufc(x) = ^(Sx), VS in the subgroup 
of SO + {n 1 1) that fixes 0. Write yu fc (|x|) to be the common value of //(x) when 

xe0B| X |(o). 

Next, we want to prove that 



/ / / sh£ _1 (x + y)dxdydv 
Js?- 1 Jo Jo 



subject to e(v) + e(— v) = E is maximum for e(v) = e(—v). This is proved using 
the augmented functional 

/ / / shf 1 (x + y)dxdy + L(e(v) + e(-v)) 
Jss;- 1 Jo Jo 

where L is the Lagrange multiplier. Setting the partial derivatives relative to 
e(v) and e(—v) to zero yields 

/ f ( ^ sh n k - 1 (e(v)+y)dy + L = 
Js^ 1 Jo 

f fe(v) 

/ / shl- 1 (x + e(-v))dx + L = 
Js?- 1 Jo 

Clearly, the above yields e(v) = e(—v) and the maximum is monotone increasing 
with E. It follows that 

But 

Mfc(0)= ,,' n „ 2 / / / sh n k -\x + y)dxdydv 

VO\(B R {0)Y Jgn-l J Q J 

Thus /Xfc(|x|) reaches its maximum at x = 0. Furthermore, it is obviously sym- 
metric for the subgroup of SO + (n, 1) that fixes 0. Finally, from the Appendix, 
/Xfc(||x||) is monotone decreasing with ||x||. 
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It remains to show that the optimal way to distribute the volume allocated 
to X is in a ball around 0. First, we show that, if mi X t (X) is achieved for 
some X* that does not enjoy the symmetry under the action of the subgroup 
of SO + (n, 1) that fixes 0, then the same Xt(X*) can be achieved for another 
subset, that has the symmetry, but that cannot be optimal. 

Decompose the big ball as Br(0) = U™L^ 1 A[r i ,r i+ i) 1 where A[ri, r i+ \) is the 
annulus {x : ri < d(0,x) < ri + \}. Clearly, there exists an annulus A[pi, P i+) C 



A [n , r i+1 ) such that A t (XnA [n , r i+1 ) ) = A t ( A[ Pi , p i+ ) ) . Thus A t (X) = £ . A t {A[ Pi , p i+ ) ) 



where T(r) is the transverse measure such that vol(A[r,r + dr)) = T(r)dr, 
that is, T(r) is the "area" of the sphere at with radius r. Consider now the 
constrained optimization problem 



Again, a Lagrange multiplier argument proves that optimality could not hold 
with some i such that Pi+ < Pi+i- (The intuition is that, since the density pu 
is monotone decreasing, optimality would require pi + i to drop down to Pi+.) 
Thus po = and p t + = Pi+\, that is, the presumed optimal annulus collapses 
to a ball. ■ 

7 Differential geometry versus real network con- 
gestion 

It is argued that, no matter how theoretical our model A t (B r (0)) of the traffic 
load is, it is remarkably accurate at predicting how the "load at the center," 
here 0, scales with the number of vertices, TV, in a real network. 

In [20], the traffic load at the "center" of a sample of networks from the 
Rocketfuel data base [22] has been numerically found to scale as ./V 2 . Using the 
scaled Gromov analysis [11], it was asserted in [20] that the networks from that 
data base are negatively curved. In [20] , the "center" is somewhat loosely defined 
as the "network core, " or the "set of nodes that are at the intersection of the 
majority of geodesies. " The latter intuitive concept is unmistakably the same 
as our theoretical center of mass concept. In the same paper [20], Narayan and 
Sanicc provide experimental evidence that the traffic load at the center of Watts- 
Strogatz Small- World 2-dimcnsional networks scales as N 1 ' 5 . It is commonly 
admitted, and it has been proved using the scaled Gromov ^-analysis [11], that 
those networks are not hyperbolic; they are rather Euclidean, even positively 



and 





subject to 
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curved, in a certain range of the connectivity and rewiring parameters [18, Sec. 
6.4.2]. 

Fact 1 (Narayan and Saniee [20]) Let (3 c (v) be the betweenness centrality 
of the vertex v in a network. 

1. In the Rocketfuel data base [22] of real networks, which are scaled- Gromov 
hyperbolic by the definition of [11], the maximum traffic rate scales as 

m&xf3 c (v) = <d(N 2 ) 

V 

2. For synthetic 2-dimensional Euclidean lattice networks, the maximum traf- 
fic rate scales as 

max/? c (u) = 0(N 15 ) 

V 

To proceed to a continuous geometry justification of the above, the graphs 
are embedded in appropriate manifolds, using the Bonk-Schramm theorem for a 
Gromov hyperbolic graph, or using the trivial Euclidean embedding for the Eu- 
clidean lattice graphs. The Bonk-Schramm embedding / :G-tflC H™ maps 
the geodesic flow on the graph G to a quasi-geodesic flow on the manifold D. 
Consider a vertex v on the graph. We clearly have A t (star(t>)) w A t (B r (f(v))), 
for some r of the order of the mesh of the graph. But as argued in Section 2, 
A t (star(w)) scales as t(v), which scales as (3 c (v). The continuous geometry model 
of the traffic metric of the above Fact is therefore A t (B r (x)) for some x G D. 
Obviously, the continuous geometry equivalent of N is 

N = vol(£>) 

Fact 2 (Theorems 2 and 3) Let A t (B r (x)) be the traffic load in a small met- 
ric ball B r (x) embedded in a large metric ball Br{Q), where R^> r. 

• In a hyperbolic ball Br(0) C H™, we have 

maxA t (B r (x)) = Q(vol(B R (0)) 2 ) = B(N 2 ) 

• In a Euclidean ball Br(0) C E n , we have 

In view of the striking consistency between Fact 1 and Fact 2, our model 
correctly predicts how the maximum traffic load scales as a function of N or the 
volume of the manifold. 

By the same token, we have extended the Euclidean results of [20] to the case 
where n > 2. This generalization allows us to observe that, as the dimension 
gets higher and higher, the Euclidean congestion decreases and the gap between 
traffic loads in hyperbolic and Euclidean spaces increases. 

While there are computational methods to check the Gromov property of 
a network [11, 18, 1, 10], associating a dimension with a complex network is 
an entirely other matter. The remarkable feature is that the asymptotic traffic 
analysis in negatively curved spaces transcends the dimension. 
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8 The case of positive curvature 



The case of a Riemannian manifold positively curved as < k\ < k < k\ is 
significantly different from that of a nonpositively curved manifold, because the 
diameter of the former is bounded as n/ki. Furthermore, by the sphere the- 
orem [7, Chap. 13], if k\ = k 2 /2 1 the manifold is homeomorphic to a sphere. 
Thus, contrary to the nonpositive curvature case, we cannot take a very large 
manifold with fixed positive curvature; so, we will have to take an arbitrarily 
large manifold of positive curvature decreasing to 0. The first part of the follow- 
ing theorem is completely trivial; the second part is easily proved by bounding 
the integrand of (4) and taking the upper limit of the two inner integrals to be 
nR = n/k. 

Theorem 6 Let M be a connected n- dimensional Riemannian manifold with 
constant curvature k = k 2 > 0. 

1. Both the normalized traffic X t and the inertia 4>^ over the whole manifold 
are uniform. 

2. X t (B r (0)) = 6 (tjjtt^), where B^ R (0) = M \ {antipodal point ofO}. 

Comparing with Theorem 3, it is clear that the normalized traffic decays with 
n even faster than in the Euclidean case. 

Finally, using the comparison formula of the Appendix, the following is easily 
derived: 

Theorem 7 Let M be a connected Riemannian manifold with curvature bounded 
as k\ < k < k\. Then 

^ 2 (x)<^^<^(x) 

Thus, under varying but bounded curvature, the traffic density in a small neigh- 
borhood of x remains bounded between the densities in fixed curvature. (Ob- 
serve that by symmetry (x), i = 1,2, are independent of x.) 



9 Conclusion 

We have provided a mathematical justification of the experimentally observed 
fact that negatively curved networks — even with uniform curvature — driven by 
uniformly distributed demand have small areas of very high traffic concentra- 
tion. Nonncgatively curved networks, on the other hand, do not exhibit this 
phenomenon as dramatically as negatively curved networks. In fact, uniformly 
positively curved networks have uniform traffic distribution; more generally, 
bounds on the curvature implies bounds on the maximum traffic. 

Since the root cause of congestion in a network is its negative curvature, load 
balancing could be achieved by controlling the curvature to become and remain 
positive [19], despite outages and varying demand. 
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The areas of maximum traffic have been narrowed down to areas of low in- 
ertia. It has been proved that networks with enough symmetry have colocated 
maximum traffic and minimum inertia. But as already said, for general net- 
works, bounding the distance between the two points is a challenging problem. 

The traffic dealt with here is the one driven by a uniformly distributed de- 
mand. The extension to nonuniformly distributed demand is not a major hurdle. 
It suffices to redefine the inertia as (/)q\v) — Y^ v i d 2 (v,v l ) (j2 j kd(v l , v^ )j . 

Finally, the present analysis is a spatial one. The temporal component would 
bring the dynamics of packet drops and retransmission into the picture. Early 
ns-2 simulation have shown that UDP traffic has its maximum packet drop at 
the point of maximum traffic/minimum inertia [19]. 



10 Appendix: Nonpositively curved spaces 

Let D be a convex domain in an n-dimcnsional Ricmannian manifold M. Here 
convex means that, for every s and t in D, there is a unique shortest geodesic, 
or shortest path, [s, t] going from s to t and contained in D. 
For a subset X of D, we are interested in 

Xt{x ^^Tnvl *(*n[s,t])dsdt 
yoL(u) j DxD 

where the integral is with respect to the square power of the n-dimensional 
volume of D. 

From a dynamical point of view, assuming uniform traffic between pairs of 
points of D, and assuming that this traffic travels at unit speed along geodesies, 
Xt(X) measures the average of the amount of time spent in X. 

We provide an estimate for Xt(X) when M has non-positive curvature, and 
an exact computation when the curvature is constant. 

Consider a point s in D. For every unit vector v based at s, draw the 
geodesic g$ emanating from s in the direction of v, and let e(v) be the distance 
from s to the point where the geodesic exits D. (This exit point is unique by 
convexity of D.) We can similarly consider the geodesic g_$ emanating from s 
in the opposite direction —v, and the corresponding distance e(—v). 

For k > 0, define 

I r re(v) re(-v) 

Mfc( s ) = yo i( D y J n _! J J p^rsinh" -1 (kx + ky) dx dy dv (3) 

where the outer integral is with respect to the (n — l)-dimensional volume of the 
unit sphere S 1 ™ -1 consisting of all unit vectors v based at s. The inner double 
integral can of course be computed by elementary calculus. 
For k = 0, let 

I r re(v) M-v) 

M s ) = im x 2 / / / (x + y) 71 - 1 dx dy dv 
VO\(DY Jgn-1 j j 
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Theorem 8 If the curvature of M is constant and equal to —k 2 , then 



X t (X) = / ^fe(x) rfx 



Jx 



More generally, if the curvature of M is everywhere bounded between two non- 
positive constants —k\ and —k\ with < fci < k^, 



In other words, the quantity Xt(X) is estimated by the integral of the density 
functions 

Proof. Let T be the set of triples (s,t,x) where s, t, x are points of D such 
that x is located on the geodesic [s,t] and inside of X. There are two natural 
measures that can be put on T. 

The first one is the product dsdtdl of the volume in s, the volume in t, 
and the arc length parameter t along the geodesic [s, t]. The reason why we are 
considering T with this measure is that 



The second one requires a different description of T. If wc consider the unit 
vector v based at x and pointing in the direction of s (so that — v points in the 
direction of t), the distance s from x to s, and the distance t from x to t, then 
the three points (s, t, x) are completely determined by the x-based vector v and 
by the numbers s and t > 0. 

Recall that the unit tangent bundle T X X consists of all unit vectors v based 
at points of X. This is a (2n — l)~dimensional manifold (n dimensions for the 
base point, n — 1 for the direction), and the metric of X naturally lifts to a 
Riemannian metric on T X X by using the Levi-Civita connection. 

The above construction identifies the set T of triples (s,t,x) with x e X n 
[s,t] to the subset of T X X x M x K consisting of those (v, s,t) where v is a 
point of T X X, and where < s < e(v) and < t < e(—v), with e(v) the 
distance between x and the point e(v) where the geodesic starting at x along 
the direction v exits D (see Fig. 3, right panel). Considering the volume form 
dv on T X X now provides another measure dv ds dt on T. 

When the curvature of M is bounded between — k\ and — fc|, the two mea- 
sures ds dt d£ and dv ds dt on T can be compared by the standard Riemannian 
arguments on the variation of geodesies (using Jacobi fields) : 

Lemma 1 





sh^" 1 (s + t) dv ds dt < ds dt dt < sh£ 2 1 (s + t) dv ds dt 
where shk(x) = ^sinh(kx) if k > and sq(x) = x. 
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Therefore, 



vol(D) 



W Jt^x Jo Jo aK;1(S + l) ^ M ^ 



f f re(v) re(-v) 

/ / / / sh™" 1 (s + t)dsdtdvdx 
J x Js™- 1 Jo Jo 



vol(D) 2 J x J s n-i J J 
Hk 2 (x) dx 



x 

Observe that, on the fourth line, dv stands for the (n— l)-dimcnsional volume on 
the sphere of unit vectors based at x, whereas dv represents the (2n — 1)- 

dimcnsional volume of T Y X in the two lines before. 
The inequality 

MX) > / Ai fcl (x)dx 
Jx 

is proved by the same argument. ■ 

For instance, consider the case where D and X are two concentric balls in 
a n-manifold of constant curvature — k 2 < 0, of respective radii R and r with 
r « R. Then, e(v) xJi for every vector v based at a point of X. 

Theorem 9 If D and X are two concentric balls of respective radii R and r, 
with r « R, in the Euclidean space of dimension n, then the proportion At(X) 
of traffic in D that transits through X is approximately 

A t (X)xc (n) ' 



where c (n) is an explicit constant depending on the dimension n. 

If, instead, D and X are concentric balls in the n-dimensional hyperbolic 
space (where the curvature is —1), then 

X t (X)^ Cl (n)r n . 

Proof. Under these hypotheses, \t(X) is of the order of /ifc(0)vol(X). Write 
/ifc(O) = Mfc(0)/vol(D) 2 , where Mfc(O) is the triple integral. e(v) x R for every 
vector v based at a point of X. Therefore, in the Euclidean case, M (0) is of 
the order of R n+1 , vol(£>) is of the order of R n , and \t(X) is of the order of 
r"/i?" _1 . In the hyperbolic case, Mx(0) is of the order of e 2 («-i)« V ol(D) is of 
the order of e ("- 1 ) i? j and thus A t (X) is of the order of vol(X), which is itself of 
the order of r" for r small. ■ 

Note the independence on R in the negatively curved case. 
Proof of Lemma: The Jacobi field equation [7] in the hyperplane orthogonal 
to the geodesic [x, e(v)] in the n-manifold M of constant curvature — k 2 reads 
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d ~ k 2 J(s) = 0. If dx,i = J(0), i = 1, ...,n— 1, is a linear clement orthogonal 
to the geodesic at x and dOi = J'(0) is the elementary angle on the sphere S 1 ™ -1 , 
the solution to the Jacobi field equation at s reads, with a similar solution for 
[x, e(— v)] at t, 

dsi = — smh(ks)d6i + cosh( ks)dxi 
k 

dti — ——sinh(kt)d9i+cosh(kt)dXi 

In the above dsi, dU are the variations of the geodesic in the orthogonal hy- 
perplanc at s, t, respectively. The volume defined by the variation of s at s 
is ds = (A™^ dsi) ds. Clearly, A™^ dsi is the transverse variation and ds is 
the longitudinal variation. A similar statement holds for dt. Because of the 
skew-symmetry of exterior differential forms, it is convenient to introduce the 
notation 6 = uj to denote 9 = {— l)™ -1 ^. With this notation, we get 

ds AdtAdi = A""] 1 (ds, A dti) Ads AdtAdi 

,n-l ' 1 



= A "=i [j- sinh(A:(s + t))d6i A cbq j Ads AdtAdi 
1 

- sinh(fc(s + tj) j {K=i de * A dx hdl Ads Adt 

= sinh"- 1 (/c(s + tj){A^d0i) A ((A^dxj) A dt) A ds A dt 

= Sinn™" 1 (k(s + t)) (dvA ((A^dxi) A dt)) A ds A dt 

= si-nh^ 1 (k(s + t))(dv A dx) Ads A dt 

= smh^ 1 (k(s + t))dv Ads A dt 

Observe that, to go from the second to the third line, we need skew-commutativity 
of A along with some elementary hyperbolic trigonometry. The inequality fol- 
lows from the Rauch comparison argument [16, VIII, Th. 4.1], which in this 
context refers to the monotone increasing property of the solution sh fe (a;) to the 
Jacobi field equation with k. ■ 

We finish this Appendix with two corollaries, useful in the main body of the 
text. A notation needs to be made more explicit: e x (v) denotes the length of 
the geodesic shot from x in the direction v and terminating at <9_Br(0). 

Corollary 1 For fixed v, e x (w) + e x (— v) is strictly monotone decreasing with 
l|x||. 

Proof. Let e x (v) be the point where the geodesic starting at x with direction 
v exits -Br(O). With a mild abuse of notation, we will let v denote the angle 
between the radial [0,x] and the tangent to the geodesic at x. Same defini- 
tion applies to e x (— v). Because HP is isotropic, the geodesic [e x (— v), e x (u)] 
is contained in a plane, which itself contains x and 0. Let [0,p] be the radial 
orthogonal to that geodesic. In the right angle triangle, A0pe x (w)), we get 
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cosh(i||e x (w)e x (-u))||) = co ^(\\ol \\ ) - 0n thc otncr nand > in thc ri S nt an S le 
triangle AxpO, we have sinh ||Op|| = cosh(x) sin v. It follows that 

, /l n , ^x.A coshi? 

cosh -||e x (w)e x (--y) 



v z / v 1 + cosh x sin v 

and from there the result is obvious. ■ 

Corollary 2 For k > 0, the density function /Ufe(s) is strictly monotone de- 
creasing with ||s||. 

Proof. It is easily seen that the double inner integral in the definition (3) of 
is a polynomial in the variables exp(±ki(e s (v) + e s (— v)), where i is some power. 
From the same definition, it is also obvious that the inner double integral is 
monotone increasing with e s (v) + e s (—v). But, by the previous corollary, the 
latter is monotone decreasing with ||s||. Thus the double inner integral in (3) is 
monotone decreasing with ||s||, and so is the average over the sphere ■ 



11 Appendix: Positively curved spaces 

The case of uniformly positive curvature k — k 2 > is treated in a way parallel 
to the preceding one. The difference resides in the Jacobi field equation, which 
now takes the form d ^Jf^ + k 2 J(s) = 0, with ordinary trigonometric function, 
rather than hyperbolic trigonometric function, solution. Therefore, if for some 
bounded subset D of the space of curvature k 2 we define 

^) = ^J^_ i ^ V) J^ ) ^^Kx + y))^ (4) 
the normalized traffic load in X C D is given by 

X t (X) = I M +(x)rfx 

In a Riemannian manifold with its curvature bounded as < k 2 < k < k%, we 
have 

/ M +(x)dx<A t P0< / /i+(x)dx 
Jx Jx 



References 

[1] F. Ariaei, M. Lou, E. Jonckeere, B. Krishnamachari, and M. Zuniga. Cur- 
vature of sensor network: clustering coefficient. EURASIP Journal on 
Wireless Communications and Networking, 2008. to appear. 

[2] Marcel Berger. Riemannian Geometry During the Second Half of the Twen- 
tieth Century, volume 17 of University Lecture Series. American Mathe- 
matical Society Providence, RI, 2000. 



21 



[3] M. Bonk and O. Schramm. Embeddings of Gromov hyperbolic spaces. 
Geom. Fund. Analysis, 10:266-306, 2000. 

[4] Martin R. Bridson and Andre Haefhger. Metric Spaces of Non-Positive 
Curvature, volume 319 of A Series of Comprehensive Surveys in Mathe- 
matics. Springer, New York, NY, 1999. 

[5] Cisco. How does load balancing work?, 2005. Document ID: 5212,2005. 

[6] M. DeVos and B. Mohar. An analogue of the Descartes-Euler formula 
for infinite graph and Higuchi's conjecture. Transactions of the American 
Mathematical Society, 359:3275-3286, 2007. 

[7] M. P. do Carmo. Riemannian Geometry. Birkhauser, Boston, Basel, Berlin, 
1992. 

[8] S. Hclgason. Differential Geometry and Symmetric Spaces. Chelsea Pub- 
lishing, Providence, RI, 2000. 

[9] Y. Higuchi. Combinatorial curvature for planar graphs. J. of Graph Theory, 
38:220-229, 2001. 

[10] E. Jonckhcere, F. Ariaci, and P. Lohsoonthorn. Upper bound on scaled 
hyperbolic gromov 5: 4-point condition. Mathematics of Networks, 1:1, 
2009. to be submitted. 

[11] E. Jonckhecrc, P. Lohsoonthorn, and F. Bonahon. Scaled Gromov hy- 
perbolic graphs. Journal of Graph Theory, 57:157-180, 2008. DOI 
10.1002/jgt.20275. 

[12] E. A. Jonckheere and P. Lohsoonthorn. A hyperbolic geometry approach to 
multi-path routing. In Proceedings of the 10th Mediterranean Conference 
on Control and Automation (MED 2002), Lisbon, Portugal, July 2002. 
FA5-1. 

[13] E. A. Jonckhcere and P. Lohsoonthorn. Geometry of network security. In 
Proceedings of the American Control Conference (ACC 2004), Boston, MA, 
June 2004. Paper WeMll.l, Special Session on Control and Estimation 
Methods in Network Security and Survivability. 

[14] J. Jost. Nonpositive Curvature: Geometric and Analytic Aspects. Lectures 
in Mathematics. Birkhauser, Basel-Boston-Bcrlin, 1997. 

[15] J. Jost. Riemannian Geometry and Geometric Analysis. Universitext. 
Springer, Berlin, Heidelberg, New York, 1998. Second Edition. 

[16] S. Kobayashi and K. Nomizu. Foundation of Differential Geometry, vol- 
ume 2. Wiley, 1996. 

[17] D. Krioukov, F. Papadopoulos, A. Vahdat, and M. Boguna. Curvature and 
temperature of complex networks. Physical Review E, 80:035101 (R), 2009. 



22 



[18] P. Lohsoonthorn. Hyperbolic Geometry of Networks. PhD 
thesis, Department of Electrical Engineering-Systems, Uni- 
versity of Southern California, 2003. Available at 
http://eudoxus.usc.edu/IW/MATTFINALTHESIS_MAIN.pdf. 

[19] M. Lou. Traffic pattern analysis in negatively curved networks. PhD thesis, 
University of Southern California, Los Angeles, CA, May 2008. Available 
at http://eudoxus.usc.edu/IW/Mingji-PhD-Thesis.pdf. 

[20] O. Narayan and I. Saniee. The large scale curvature of networks. 
arXiv:0907.1478vl, [cond-mat.stat-mcch] 9 July 2009, July 2009. 

[21] T. Rcti, E. Bitay, and Z. Kosztolanyi. On the polyhedral graphs with 
positive curvature. Acta Polytechnica Hungarica, 2(2):19-37, 2005. 

[22] N. Spring, R. Mahajan, D. Wetherall, and T. Anderson. Measuring ISP 
topologies with rocketfuel. IEEE Transactions on Networking, 12(1):2-16, 
February 2004. 

[23] Liang Sun and Xingxing Yu. Positively curved cubic plane graphs arc finite. 
Journal of Graph Theory, 47(4):241-274, 2004. 

[24] L. Zhao, Y.-C. Lai, K. Park, and N. Ye. Onset of traffic congestion in 
complex networks. Physical Review E, 71:026125-1-026125-8, 2005. 



23 



